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The hard X-ray view of the low-luminosity blazar 
in the radio galaxy NGC 6251 
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Abstract. We present results from a BeppoSAX (July 2001) observation of the FR I radio galaxy NGC 6251, 
together with a re-analysis of archival ASCA (October 1994) and Chandra (September 2000) data. The weak 
detection above 10 keV and the lack of iron fluorescent emission lines in the BeppoSAX spectrum rule out 
that the bulk of the X-ray emission is due to an obscured Seyfert nucleus. The study of the multiwavelength 
spectral energy distribution suggests instead that X-rays probably originate as inverse-Compton of synchrotron 
seed photons in a relativistic jet, indicating that NGC 6251 hosts a low radio luminosity (L5 ghz ~ 10*" erg s~^) 
blazar. The BeppoSAX spectrum is flatter than in the earlier ASCA observation. This might be due to the 
emergence of a different spectral component during phases of lower X-ray flux. In this context, we discuss some 
possible explanations for the intense and mildly-ionized fluorescent iron line measured by ASCA. 
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1. Introduction 

; NGC 6251 is an E2 nearby (z = 0.02488; if Ho = 
. 70 km s~^ Mpc~^, as assumed in this paper, 1" 
' corresponds to about 500 pc) radio galaxy with a 
Faranoff-Riley I morphology, known to host a gi- 
^ ant r adio jet (|Wagget et al. 1977| ), a Seyfert 2 nu- 
cleus IjWerner et al. 2002|l . and almost edge-on dust lanes 
l|Nieto et al. 1983|l . The nucleus is likely to contain a su- 
permassive black hole with mass M 4-8 x IQ^Mq 
HFerrarese fc Ford 1999|l . as suggested by the presence of 
nuclear gas and a dust disc on scales of a few hundreds 
parsecs. NGC 6251 belongs to the outskirts of the cluster 
Zw 1609.0-F8212 ( |Young et al. 1979| ), whose influence on 
the properties of the galaxy should be, however, marginal 
UPrestage fc Peacock 1988) 1. 

The NGC 6251 jet is one of the most spectacular radio 
objects of the whole sky. It is a 4'. 5 long, highly coUi- 
mated (opening angle 7.4°) structure ) |Perley eTal. 1984) . 
Regions of enhanced radio emission along the jet were la- 
beled by Birkinshaw & Worrall (1993) in the 330 MHz 
radio map as A (10-40" from the nucleus), B (40-126"; 
actually structured in smaller sub-structures) and D (178- 
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264", where the jet bends towards the North; it is as 
well highly structured). Between B and D one finds a re- 
gion of low radio brightness (C). None of these region 
was observed to be a significant source of X-ray emis- 
sion by ROSAT (Birkinsh aw fc Worrall T993|l . A sub-pc 
counter-jet was discovered only recently with Very Long 
Baseline Interferomer observations (jSudou et al. 200 
U-band Hubble Space Telescope images unveiled a region 
of extended emission, lying nearly perpendicular to the 
radio jet axis and the dust ring IjCrane fc Vernet 1997|l . 
which is likely to originate from scattering of a nuclear 
continuum source. 

Discovered in the X-rays by the Einstein 
Imaging Proportional Counter IjJones et al. 1986)l . 
NGC 6251 showed in the RO SAT PSPC observa- 
tion l|Birkinshaw fc Worrall 1993|l an unresolved core 
{FWHM < 4"), with a possible extended halo on scales 
<; 100 kpc IjMack et al. 1997,) . The first observation 
of NGC 6251 covering the intermediate X-rays (i.e.: 
2-10 keV) was performed by ASCA. The presence of a 
bright Kq, fluorescence iron from ionized iron (centroid 
energy, - 6.68 keV; Equivalent Width, EW ~ 600 eV; 
ITurner et al. 19971 ISambruna et al. 1999)l . and the fact 
that the continuum could be best fit with a combination 
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of a standard AGN power-law component (photon index, 
r, of 2.11) plus a thermal soft excess, suggested that ion- 
ized gas may significantly contribute in this energy band. 
Electron scattering of an otherwise invisible nuclear con- 
tinuum was an interesting possibility, in light of the HST 
discovery of extended ionized gas, which could potentially 
act as a nuclear mirror ("warm mirror" hereinafter). 
Recently, an association has been proposed between 
NGC 6251 and the EGRET source 1EGJ162H-8203 
HMukherjee et al. 2002| ). If this is confirmed, NGC 6251 
would be one of the three radio galaxies - alongside 
with Cen A (Srecku^ ^et al. 199911 and 3EG J1735-1500 
IjCombi et al. 2003) - detected in high-energy 7-rays. No 
detection by the Extreme Ultraviolet Explorer is reported 
in the literature UMarshall et al. 1995|l . 

The deep BeppoSAX observation, described in this 
paper, aimed at verifying the interpretation scenario 
emerging from the ASCA outcomes. Thanks to its un- 
surpassed sensitivity in hard X-rays (i.e.: >10 keV; 
IBoella et al. iggTji . BeppoSAX would be easily able to de- 
tect a transmitted nuclear component piercing through a 
Compton-thick absorber IjMatt et al. 2000|l . The results of 
this observation are presented in this paper, together with 
an analysis of archival ASCA and Chandra data of the 
same source. The log of the observations discussed in this 
paper is presented in Table 

While this paper was in the process of being ac- 
cepted, we became aware of a recently accepted paper by 
Chiaberge et al. (2003), who discuss the same X-ray ob- 
servations, together with literature and HST data. Their 
X-ray results are largely coincident with ours. 

2. BeppoSAX results 

2.1. Data reduction 

BeppoSAX data were reduced according to standard pro- 
cedures as in, e.g., Guainazzi et al. (1999). Scientific 
products for the imaging Low Energy Concentrator 
Spectrometer (LEGS, Parmar et al. 1997. 0.5-4 keV) 
and Medium Energy Concentrator Spectrometer (MEGS, 
IRoella et al. 1997bl 1.8-10.5 keV) were extracted from cir- 
cular regions of 8' and 4', respectively. Background spectra 
were extracted from blank sky field event lists provided by 
the BeppoSAX Science Data Center, and appropriate for 
the date of the observation. Products for the Phoswitch 
Detector System fPDS. ITVoirtera et al. 19971 13-200 keV) 
were extracted by plain subtraction of those corresponding 
to intervals when the instrument was observing NGC 6251, 
and a region ±3.5° aside, according to the standard 96 s 
cycle. 

All the spectra employed in this paper have been re- 
binned in order to oversample the intrinsic energy reso- 
lution of the instruments by a factor not larger than 3, 
and to have a number of counts in each spectral channel 
higher than 30, in order to ensure the applicability of the 

test. In this paper: energies are quoted in the source 
rest frame; uncertainties on the spectral parameters are 
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Fig. 1. BeppoSAX spectrum {upper panel) and residuals 
in units of standard deviations {lower panel), when the a 
photoelectrically absorbed power-law model is applied. 



quoted at the 90% confidence level for one interesting pa- 
rameter; uncertainties on the count rates are at Icr level, 
unless otherwise specified. 

2.2. Spectral results 

In Fig.Hthe results of the fit of the NGC 6251 BeppoSAX 
spectrum with a simple absorbed power-law is shown. 
The fit is statistically acceptable (x^ = 96.7/92 degrees 
of freedom, dof). The best-fit parameters and results 
are reported in Table [3 The PDS data points lay ac- 
tually above the best-fit model, which is mainly deter- 
mined by the statistics of the imaging instruments. The 
excess is, however, not highly significant. The ratio be- 
tween the PDS flux and the extrapolation of the best 
fit model in the LEGS/MEGS energy band is 1.4 ± 0.6, 
therefore only marginally exceeding the typical values 
of the cross-normalization factors between the PDS and 
the MEGS f0.80-0.85: lFiore et al. 1998}! . Not surprisingly, 
several models provide comparably good descriptions of 
this hard excess, but none of them is required from the 
statistical point of view. A broken power-law, for instance, 
yields a Ax^ = 3.4 for a reduction of the number of degrees 
of freedom {Av) by 2 (this quantity will be indicated as 
Ax^ /Aiy hereinafter), corresponding to a confidence level 
-83.0% 

The BeppoSAX large source spectrum extraction re- 
gion encompasses the putative ROSAT extended emis- 
sion. We therefore tried to add a soft X-ray compo- 
nent, modeled by a coUisionally ionized plasma (we used 
the mekal implementation in XsPEC Vll.O through- 
out this paper; I Mewc et al. 19851 ILiedahl et al. 19951 
[Arnaud fc Rayniond 1992) . Again, the improvement in 
the quality of the fit is not significant {Ax'^/Aiy — 4.4/2). 
Consequently, the temperature is rather poorly deter- 
mined {kT = I.4±0:f). The 0.5-2.4 keV flux of this com- 
ponent would be (1.3 ± 0.7) x 10^^^ erg cm^^ s^^, cor- 
responding to about 15% of the non-thermal component 
in the same energy band. This model requires an absorb- 
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Satelhte 


Start time 


Exposure time 


Count rate 






(ks) 


(10^^ s"^) 


ASCA 


28-Oct-1994 


36° /44'' 


7.2 ±0.2"/3.63± 0.13" 


Chandra/ AClS-l 


ll-Sep-2000 


25 


8.1 ±0.2 


BeppoSAX 


19-Jul-2001 


28780'*/72'= 


2.78 ± 0.1275.55 ± 0.09''/13 ± 4= 



"SISO, ''GIS2, '^LECS, ''MEGS, "^PDS. 

Table 1. Log of the observations presented in this paper. 
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Nh 
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fcTi/Fsoft 


fcT2 /inbreak 
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EW 


XV dof 


Mission 


(10^^ cm-2) 




(keV)/ 


(keV) 


(keV) 


{IQ-^ ph cm-2 s-i) 


(eV) 




BeppoSAX 


















PL 


o.8±g:i 


i.79±g:g? 












96.7/92 


PL+TH 


i.o±g:Si 


i.70±E!:i^ 


i.4±g:? 




6.4I" 
6.59I" 
6.7I" 
6.961" 


< 8.0 

< 4.5 

< 3.0 

< 2.1 


< 160 

< 80 

< 60 

< 50 


92.3/90 


ASCA 


















PL 


o.9±g:| 


2.06±!];E!? 












259.7/225 


PL+TH 




2.oo±g:gf 


o.85±8:f| 










240.6/223 


PL+2xTH 


1.6±gi 


2.3 ±0.4 


o.8±[!:i 


6±1 








231.8/221 


PL+TH+GA 


J-- ' ±0.4 


2.05±S]:|,J 






6.59 ±0.16 


8±4 


500 ± 200 


228.1/221 


BKPL+GA 


4.o±l:i 


2.28±o:i4 


3.7±;:i 


1.38+0:^9 


6.62±E!:li 


10 ±4 


800 ± 300 


245.0/226 



1^ Gaussian emission profiles added to the best-fit continuum. Centroid energies are fixed. Upper hmits are at the 90% 
confidence level 

Table 2. Best-fit parameters and results for the analysis of the nuclear spectra of NGC 6251. Model legenda: 
PL = power-law; TH = thermal component (mekal implementation in Xspec) with solar abundances; GA = Gaussian 
emission line; BKPL = broken power-law. 



ing column density slightly in excess with respect to the 
Galactic contribution along the line-of-sight to NGC 6251 
(5.7 X 10^° crn' MMlIT^hy et al. 1996| ). 

No systematic residual feature is present at the en- 
ergies, where Kq fluorescence transitions of iron are ex- 
pected. Negligible improvements of the are yielded by 
the addition of a narrow (i.e., intrinsic width, a — Q) 
Gaussian profile to the power-law model. 90% upper limits 
on the intensity of a neutral (6.4 keV) or He-like (6.7 keV) 
iron lines are 8.5 and 3.0 x 10^® photons cm^^ s^^, re- 
spectively, corresponding to EWs of 160 and 60 eV, re- 
spectively. These upper limits are largely inconsistent with 
the ASCA detection HTurner et al. 19971 cf. also Sect. 3.1 
later). Simulations show that a 600 eV EW line would have 
been detected at the 6-7cr confidence level in the MEGS 
spectrum. 

The 0.5-10 keV flux during the BeppoSAX observa- 
tion was 7.5 X 10~^^ erg cm~^ s""'^, corresponding to un 
unabsorbed rest-frame luminosity of 1.13 x 10"*^ erg s~^ 
in the same energy band. 

3. Comparison with previous X-ray observations 

3.1. ASCA 

We have retrieved the data of the ASCA observation of 
NGC 6251 from the public archive as screened event lists. 



Spectra were extracted from regions of radii 4', 3'.1 and 
3'. 75 in the SISO, SISl (grade 0, 2, 3, 4) and GIS detectors, 
respectively (they therefore encompass as well the ROSAT 
extended emission region). Background spectra were ex- 
tracted from regions in the field of view of the detectors, 
free from contaminating sources. Response matrices were 
generated with the packages included in the Lheasoft 
v5.0 software. The spectral analysis was performed in the 
0.5-10 keV and 0.7-10 keV energy bands, for the SIS and 
GIS instruments, respectively. The results of our analysis 
substantially coincide with those presented by Turner et 
al. (1997) and Sambruna et al. (1999), and we summarize 
them in this paper for the sake of clarity only. 

A simple power-law is a fairly good representation of 
the ASCA spectra (x^ = 259.7/225 dof). However, the 
addition of a mekal component improves significantly the 
quality of the fit (Ax^/Ai/ — 19.1/2, corresponding to 
an F-test confidence level ci99.97%). The addition of a 
multitemperature blackbody (model diskbb in Xspec) 
yields, on the contrary, a negligible improvement to the 
quality of the fit. Modeling the continuum in terms of 
a broken power-law yields a comparatively worse fit as 
well. A systematic excess around the energy where fluo- 
rescent Kq, transitions from iron are expected is observed 
(see Fig. El). It can be cured with an unresolved Gaussian 
profile (Ax^/Az/ = 12.5/2, corresponding to a confidence 
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Energy (keV) 

Fig. 2. GIS2 {filled dots) and GIS3 {open circles) residuals 
against a power-law plus optically thin thermal plasma 
continuum in the 4-10 keV bands. 

level ~99.7%), or with a hotter {kT ~ 6 keV) thermal 
component (Ax^/Az/ = 8.8/2, corresponding to a confi- 
dence level ~98.6%). 

In the best-fit model (cf. Table |2Jl, the power-law in- 
dex (r ~ 2.06) is significantly steeper than observed by 
BeppoSAX. The temperature of the colder thermal com- 
ponent is kT ~ 0.8 keV, and its 0.5-2.4 keV flux [(1.6 ± 
0.4) X 10""'^'^ erg cm~^ s~^] is in good agreement with that 
of the corresponding component in the BeppoSAX spec- 
trum. The iron line centroid energy is consistent with Kq, 
fluorescence from mildly ionized to He-like iron. Its EW 
is 500 ± 200 eV, significantly larger than the upper limits 
measured by BeppoSAX. The difference is most likely due 
to a change in the underlying continuum, as the intensity 
of the ASCA line [(8 ± 4) x 10"^ photons cm^^ s^i] is 
marginally consistent with the BeppoSAX upper limits. 
However, an intrinsic variability of the iron line cannot be 
ruled out. 

The observed 0.5-10 keV flux during the ASCA obser- 
vation was 2.7 X 10~^^ erg cm~^ s~^, corresponding to an 
unabsorbed total luminosity of 4.4 x 10'*^ erg s^^ in the 
same energy band. 

3.2. Chandra 

Data of an ACIS-I NGC 6251 Chandra observation were 
retrieved from the public archive. A bright source [count 
rate (8.1 ± 0.2) x 10~^ s~^] is detected with coordinates: 
0^2000 = 16''32™3P.8; (Jaooo = +82''32™16^ i.e. 0.6" dis- 
tant from the optical nucleus of NGC 6251, well within 
the accuracy of the Chandra attitude reconstruction. At 
this count rate level, a Chandra source is likely to be 
substantially affected by pile-up, given the instrumental 
mode employed (Time Exposure Mode with a 3.2 s frame 
time). We have tried to fit the 0.5-8 kcV spectrum ex- 
tracted from the innermost 5" with a photoelectrically 
absorbed power-law, corrected for pile-up according to 
the XSPEC implementation of J.Davis' algorithm (model 



J. l_,CtilvJil 


d 


Area 


CR 


p 






(arcsec) 


(arcsec'^) 


(10-*) 


a 


(nJy) 


A 


10-40 


46.8 


2.8 


2.3 


0.34 


B 


40-126 


140.6 


2.2 


1.8 


0.26 


C 


126-178 


78.1 


2.9 


2.4 


0.36 


D 


178-264 


124.9 


1.9 


1.6 


0.24 


the 0.5-10 keV band, 


in units of 10 er; 


5 cm" 


-2 s-i 



Table 3. 90% upper limits on the ACIS-I count rates 
{CR), 0.5-10 fluxes (F), and 1 keV flux density {S^) 
across the NGC 6251 jet regions (following their defini- 
tion in Birkinshaw & Worrall 1993). The distance d is 
defined from the inner-outer border of each region to the 
core. Fluxes are calculated assuming a power-law spec- 
trum, with r = 2 and photoelectric absorption column 
density Nh = 5.7 x 10^° cm'^. 



# 


a20oo 


(52000 


count rate 
(10-* s-i) 


1 


16''32'"31'.8 


+82''32'"16'' 


810 ± 20" 


2 


16''33™00^5 


+82''31'"12'' 


4.8 ± 1.8 


3 


16''33™21^3 


+82''31'"56'' 


5.6 ± 1.8 



"NGC 6251 nucleus 

Table 4. Sources detected in the ACIS-I observation of 
the NGC 6251 field, within 4' from the core. 



pileup in XsPEC; lArnaud &: Dorman 2002|l . The fit is 

acceptable (x^ = 57.1/78 dof), and yields best-fit pa- 
rameters which are consistent with those measured dur- 
ing the BeppoSAX observation, save a 60% lower flux: 
Nh = (1.6±0.5) X 10^1 cm-^; T = 1.76±0.16; 0.5-10 keV 
flux of 4.5 X IQ-^^ erg cm-^ s-^ 

Thanks to the unprecedented spatial resolution of the 
ACIS-S, the issue of the spatial extension can be bet- 
ter addressed. No evidence for extended emission along 
the arcminutes scale jet is detected, as the ROSAT ob- 
servation had already shown in the soft X-rays only 
UBirkinshaw fc Worrall 199311 . The upper limits on the X- 
ray fluxes of the regions A to D are reported in Table 01 
They are by a factor 5 to 10 tighter than those determined 
by Birkinshaw & Worrall (1993). However, the upper limit 
on the knot "D" flux density is inconsistent with the de- 
tection {Su = 13 ± 2 nJy) obtained by Mack et al. (1997) 
from a reanalysis of the same ROSAT/PSPC observation 
discussed by Birkinshaw & Worrall (1993). 

In the innermost 4' around the NGC 6251 core only 
two sources, alongside with the nucleus itself, are detected 
at a signal-to-noise ratio higher than 3. Their positions 
and count rates are reported in Table 01 None of them 
exhibits a clear association with any radio structures. 
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Fig. 3. PDS count rates in the 13-200 keV energy band 
against the column density of a Compton-thick absorber 
covering the NGC 6251 nucleus, for different values of 
the scattering fraction /. The shaded area indicates the 
PDS detection yielded by the BeppoSAX observation of 
NGC 6251. / is defined by the ratio between the nor- 
malizations of the transmitted and the warm scattered 
components, assuming an optically thin scatterer 
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Fig. 4. Radio (6 cm). X-ray (0.5-4.5 keV) and IR (25^) 
color-color diagram. The dotted line represents the blazar 
track according to the blazar unified scheme of Fossati et 
al. (1998). The filled circles represent the data points cor- 
responding to the NGC 6251 ASCA and BeppoSAX ob- 
servations. The cross corresponds to the Seyfert 2 sample 
of Mass-Hesse et al. (1995). The empty square represents 
the SED of radio-quiet quasars IjElvis et al. 1994|l 



4. Discussion 

4.1. The nature of the nuclear X-ray emission in 
NGC 6251 

The main goal of the BeppoSAX observation described in 
this paper was to test whether the strong (EW ~ 600 cV), 
ionized Kq, fluorescent iron line observed in the ASCA 
spectrum HTurner et al. 1997|l was due to a "warm mir- 
ror" reflection-dominated, Compton-thick Seyfert 2 spec- 
trum. In the light of the BeppoSAX observation out- 
comes, this possibility is rather unlikely. The PDS de- 
tection (0.13 ± 0.04 s-i in the 13-200 keV energy band) 
is statistically consistent with the extrapolation of the 
2-10 keV spectrum. This constraints the column den- 
sity of any Compton-thick absorber to the nucleus to be 
^ 3 X 10^4 (6 X 10^4) cm-2 for a 10% (1%) scattering frac- 
tion (see Fig.OJ. Moreover, the 2-10 keV BeppoSAX spec- 
trum lacks any trace of fluorescence ionized iron lines, 
which are almost always observed in reflection-dominated 
Seyfert 2 galaxies HTurner et al. 1997bllM'att et al. 2000|l . 
The upper limits on their EWs are rather strict, ranging 
between 60 and 80 eV for He- and H-like iron. The stan- 
dard AGN 2-10 keV spectral index (F ~ 1.7), and the 
lack of Ka neutral iron line (again, with a rather tight 
EW upper hmit: 160 cV) in the BeppoSAX spectrum rule 
out a significant contribution by Compton-reflection from 
the inner side of the molecular torus ("Kr olik et al. 19941 
IGhisellini et al. 19941 IMatt et al.~19961 In principle, the 



line flux could be suppressed by resonant trapping at the 
ionized surface of an accretion disk IjMatt et al. 19931 . 
However, further pieces of evidence rule out that the X- 
ray spectrum is dominated by reflection of an otherwise 
invisible Seyfert nucleus. 

Some hints come from the study of the multiwave- 
length Spectral Energy Distribution (SED). Recently, 
Fossati et al. (1998) have proposed a unified scheme to 
explain the multivawelenght SED of blazars, whose prop- 
erties (peak frequency of the synchrotron and inverse- 
Compton components, luminosity ratio between them) are 
mainly governed by a single parameter related to the over- 
all luminosity. This scheme can be applied to help iden- 
tifying the nature of the bulk of the X-ray emission from 
the NGC 6251 nucleus. In Fig. 0]we show a color-color 
diagram between the radio (ly = 6 cm), soft X-ray (0.5- 
4.5 keV), and infrared (A = 25/i) luminosities. The dot- 
ted line in this diagram represents the trace followed by 
blazars according to the Fossati scheme. In the same plot 
the cross represents the position of the Seyfert 2 galax- 
ies of the Mass-Hesse et al. sample (1995). The 6 cm ra- 
dio luminosity of N GC 6251 is (1.0 ± 0.1) x 10*° erg s^^ 
Ij Jones et al. 1986)l . The data points corresponding to the 
ASCA and BeppoSAX observations lay intriguingly well 
on the blazar track. This supports the idea that the bulk 
of the nuclear emission in NGC 6251 is due to a jet. 

Chiaberge et al. (1999) discovered a clear correla- 
tion between the optical (Fq) and the radio core lumi- 
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nosity (Ffj) in a sample of FRI galaxies extracted from 
the 3C catalogue and observed with the HST WFC2. 
Assuming the NGC 6251 SED published by Ho (1999), 
\og{Fji/Fo) — 3.4, in perfect agreement with the value 
derived from the 3C correlation (3.7 ± 0.4; the uncertain- 
ties represent the r.m.s. scattering of the data points in the 
correlation), and largely inconsistent with values typically 
measured in radio-loud quasars (—-1; lElvis et al. 1994^ . 
This supports a common origin for the radio and optical 
emission as synchrotron radiation. HST observations in 
the U-band measured a rather high degree of polarization 
(close to 50%) in clumps close to the nucleus along the 
radio axis. Crane & Vernet (1997) suggest that the UV 
emission of these clumps is due to scattering. However, 
such an evidence is consistent as well with the possibility 
that the synchrotron-dominance extends well within the 
UV range. 

We tried to fit the overaU NGC 6251 SED from radio to 
7-rays (X-rays represented by the BeppoSAX spectrum) 
with an homogeneous Synchrotron Self-Compton (SSC) 
model l|Tavecchio et al. 1998|l . This model assumes that 
synchrotron radiation is produced by relativistic electrons 
with density n,, continuously injected in a spherical region 
of radius r with a magnetic field B and moving with bulk 
factor Fbuik at an angle 9 with respect to the line of sight 
(the Doppler factor is therefore S — [Fbuik(l— /3cos0)]^^). 
These photons are subsequently upscattered by the same 
electrons. Following Tavecchio et al. (1998), the electron 
energy distribution is modeled with a broken power-law 
with indices ni < 3 and n2 > 3 below and above a break 
energy 7^, respectively. The IRAS data points, although 
in principle available, have not been included in the fit 
due to the unknown contamination from the host galaxy. 
The best fit from radio to 7 frequencies (dashed line in 
Fig. EJ fails to reproduce the BeppoSAX spectral shape 
and slightly deviates below 10^° Hz. If one forces the 
model to match the BeppoSAX data point [dotted line 
in Fig.|SJ), the radio emission is more overproduced. This 
may indicate that self-absorption (which is not explicitly 
included in SSC models) plays an important role redwards 
the synchrotron peak. Given the fact that the SED data 
points correspond to non-simultaneous observations, one 
cannot push too far the comparison between the model 
and the observation. Nonetheless, the overall properties of 
the best-fit SSC models are not strongly dependent on the 
details of the fit. The synchrotron peak is below 10^'* Hz, 
whereas the Compton peak remains constrained between 
lO^"'^'^ and 10^^'^ Hz. As a reference the best-fit parameters 
for the fits are in the range: B = 0.01-0.15 G, S = 2-4, 
m = 1.75-2.6, n2 = 4.1-5.4, 7^ = 6.3 x 10^-1.6 x lO'' 
rig = 10^-10^ cm~^, and r = 1-8x10^^ cm. The reader 
is referred to Chiaberge et al. (2003) for a more detailed 
comparison between the NGC 6251 SED and SSC models. 

It is hard to interpret the spectral variability between 
the ASCA (empty circles in Fig. (S)) and the BeppoSAX 
spectra in terms of pure SSC model. In principle, the 
ASCA spectrum may represent the trailing edge of the 
synchrotron component. However, this would imply a shift 
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Fig. 5. NGC 6251 SED. Data are from not-simultaneous 
observations compiled by Ho (1999), save the and X- 
ray data (BeppoSAX, dots; this paper), and the putative 
EGRET detection. The dashed line represents the best 
radio-7 fit SSC model; the dotted line represents the fit 
with an SSC model, forced to reproduce the BeppoSAX 
spectral shape. The ASCA spectrum (empty circles) is 
shown as well for reference 

of the synchrotron peak by more than 4 orders of magni- 
tude, which is rather unlikely. Alternatively, the steeper 
ASCA spectrum may imply that a different spectral com- 
ponent may be emerging during phases of low X-ray 
flux. This component may be completely outshined dur- 
ing BeppoSAX- like. X-ray brighter states. We will further 
discuss this possibility in Sect. 4. 2, in connection with the 
strong fluorescent iron line observed in the ASCA spectra 
only. 

Recently, it has been suggested that the bulk of the 
X-ray emission in radio galaxies may be due to a hot 
inner accretion flow, following a line of thought sug- 
gested more than 20 years ago by Rees et al. (1982). Ho 
(1999) estimates the bolometric luminosity of NGC 6251 
as ~10^^ LEdd- This may indicate that an Advection 
Dominated Accretion Flow (ADAF; [Narayan fc Yi 1995| ) 
is responsible for the bulk of the emission in the 
NGC 6251 core. In this scenario. X-rays are mainly 
produced via bremsstrahlung by a distribution of ther- 
mal electrons with typical temperatures kT ^ 100 keV 
HDi Matteo et al. 2 000). In Fig. we compare the ratio 
between the V and the 1 keV flux density (BeppoSAX 
measurement) in NGC6251 with the values expected by 
the blazar unification scenario HFossati et al. 1998|l and 
by models of ADAF applied to the multiwavelength SED 
of nearby elliptical galaxies suspected to host low radia- 
tive efficiency accretion flows IjDi Matteo et al. 20001 we 
consider hereby models with wind, and outermost hot ac- 
cretion radius of 300 gravitational radii). The NGC 6251 
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Fig. 6. Radio between the 1 keV and the V-band flux 
density in: a) NGC 6251 {filled circle); b) the blazar radio 
luminosity classes in the unified scenario after Fossati et 
al. (1998); c) ADAF models applied to a sample of nearby 
elliptical galaxies (|Di Matteo et al. 20001 details in text). 
In the last case, the error bars on the y-axis represent 
the r.m.s. of the sample values in the corresponding radio 
luminosity interval. 



data point is in principle consistent with both scenarios. 
Applying a bremsstrahlung model to the BeppoSAX spec- 
trum, one gets indeed a total X-ray luminosity Lx / L^dd — 
(2 ± 1) X 10"'^. However, the fit is significantly worse 
than with a simple power-law (x^ = 110.0/92 dof). A 
composite bremsstrahlung and power-law fit (the latter 
component taking into account possible Comptonization 
effects) yields an implausibly low electron temperature 
{kT ~ 3 keV). These pieces of evidence are in agreement 
with the conclusions drawn by Ferrarese & Ford (1999), 
who remark that the nuclear non-thermal bolometric op- 
tical luminosity in NGC 6251 is larger than predicted for 
accretion at the Bondi rate with 10% efficiency, making 
the case for an ADAF far less compelling. 

4.2. Additional components in the nuclear X-ray 
spectrum of NGC 6251 

As already noticed by Turner et al. (1997), a soft X- 
ray excess is observed in the NGC 6251 ASCA spec- 
tra. A multitemperature blackbody is not able to fit the 
excess, arguing against the possibility that this compo- 
nent originates in an accretion disk. The best-fit tem- 
perature of the soft excess, if modeled with an opti- 
cally thin, coUisionally excited plasma, is kT ~ 0.8 keV, 
with an unabsorbed luminosity of ~ 3.6 x lO''^ erg s^^. 
Such a temperature is typical of gaseous halos in ellipti- 
cal galaxies ()Matsumoto et al. 1997|l . The parameters of 



this component are consistent with those of a similar 
component in the BeppoSAX spectrum, which is, how- 
ever, not required from the statistical point of view. The 
ASCA temperature is in turn consistent, within the rather 
large statistical uncertainties, with the temperature of 
the extended halos measured by the Einstein/WG and 
the ROSAT/ PSPC in NGC 6251 (for th e latter, kT = 
0.3-0.8 keV; IBirkinshaw <k Worrall 19*931 and references 
therein). Similarly, the 0.5-2.4 keV unabsorbed fluxes 
of this component measured by ASCA and BeppoSAX 
are consistent with the ROSAT/PSPC soft X-ray "halo" 
flux, integrated across its whole extension [(2.4±2;§) x 
10-13 erg cm-2 3"^ IMack et al. 1997j . It is therefore 
straightforward to identify the soft excess in the large 
ASCA aperture with the diffuse emission in ROSAT. Mack 
et al. (1997) discuss (and rule out) the possibility that 
such a plasma could be the conflning medium for the jet, 
along the whole structure up to knot "D" . A ~10^ K ther- 
mal emission responsible for the confinement of the jet in 
NGC 6251 should not extend beyond scales larger than 
60 pc. On such small scales, it might be marginally resolv- 
able by Chandra. Unfortunately, the pile-up affecting the 
Chandra observation prevented us from deriving precious 
constraints on the X-ray extension around the NGC 6251 
core. 

The ASCA spectrum unveils the possible presence of 
a third spectral component, whose main signature is a 
bright {EW ~ 600) eV and significantly ionized Kq, flu- 
orescence iron line. This component, even if present with 
comparable flux in the BeppoSAX observation, would 
have easily missed detection, due to the fact that the 
BeppoSAX non-thermal continuum 6 keV flux density 
was >6 times larger than in ASCA. This "third" com- 
ponent cannot be produced by the superposition of un- 
resolved discrete sources, integrated in the large ASCA 
aperture. The Chandra image shows only two sources in 
the innermost 4' around the NGC 6251 core, whose to- 
tal count rate [(1.1 ± 0.3) x 10^'^], corresponds to a 0.5- 
9 keV flux ~1.0 X 10"^"* erg cm~^ s~^ (assuming a thermal 
plasma with kT = 6 keV). This is more than one order 
of magnitude less than the flux of the ASCA component 
in the same energy band (2 x 10""'^'^ erg cm"'^ s~^). This 
component may therefore represent the "tip of the ice- 
berg" of an underlying Seyfert nuclear "warm scattering" 
(Matt et al. 1996, 2001), whose relative contribution be- 
comes not- negligible during phases of lower blazar activity. 
Alternatively, this component may represent a hotter and 
weaker phase of the thermal emission. It has already been 
pointed out that the kpc-scale jet structure in NGC 6251 
requires an ambient gas with temperatures of 2-5 keV to 
conflne it. So far, the lack of detection of a gas component 
with such an high-temperature had led to the problem 
of an "over-pressurized" jet in NGC 6251. It would be 
tempting to speculate that the high-temperature thermal 
component measured by ASCA is the required "missing 
link" to ensure the thermodynamical stability of the kpc- 
scale jet. Unfortunately, only rather strict upper limits can 
be set on any X-ray emission associated with the radio jet 
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(see Table 13). Alternatively, the line could be produced in 
the interaction between the "bloated" base of a jet with a 
stratified velocity structure and the circumnuclear matter, 
as proposed by Chiaberge et al. (2000) in the context of 
a possible unification model between BL Lac objects and 
FR I radio galaxies. The possibility that the iron line is 
associated with the highly ionized "skin" of an accretion 
disk cannot be in principle ruled out, although evidence 
for significantly ionized disk lines is rather scanty so far. 

All the best-fit X-ray models discussed in this pa- 
per require an amount of cold photoelectric absorption 
in excess to the Galactic contribution. The weighted be- 
tween the BeppoSAX and ASCA measurements is Nh = 
(1.06 ± 0.11) X 10-^^ cm^^. This is in excellent agree- 
ment with the column density through the dusty disk, 
as derived from its visual extinction {Ay = 0.61 ± 0.12; 
IFerrarese &: Ford 1999|l . if standard gas-to-dust Galactic 
ratios are assumed. This is in agreement with the idea that 
the standard pc-scale optically and geometrically thick 
torus - even if present in FRI low-luminosity radio galax- 
ies - does not intercept our line-of-sight to the nucleus 
UGhiaberge et al. 2002| ). 
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